p=struct("rho",0.8,"Ts",1.5,"eps",0.1);
time=0:0.01:20;
mus=time./((1-(p.rho)*time+p.rho*p.Ts)*p.eps);

figure;hold on;grid on
plot(time,mus);
plot([min(time) max(time)],1/p.eps*ones(1,2));
ylim([-20 20])